HIGH SENSITIVITY ANTI-NEUTRINO DETECTION BY KAMLAND 
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Electron anti-neutrinos (i7 e ) from nuclear power reactors have been observed by KamLAND. 
Data from 0.16 kton-year exposure (145.1 live days) indicates disappearance of i7 e at 99.95% 
C.L. in the energy range 2.6 MeV < £„- < 8.0 MeV. Considering two-flavor neutrino oscillation 
with CPT invariance, the only remaining solution to the solar neutrino problem is the Large 
Mixing Angle (LMA) solution. In addition a 0.28kton-year exposure (185.5 live days) was 
searched for i7 e in the energy range 8.3 MeV < E,r e < 14.8 MeV. No candidate events were 
found with expected background of 1.1±0.4 events. Assuming that the origin of v e in this 
energy region comes from 8 B solar v e , we find an upper limit of v e to v e conversion probability 
of 2.8 xl0~ 4 . 



1 KamLAND Experiment 

1.1 Motivation of the Experiment 

Kamioka Liquid scintillator Anti-Neutrino Detector (KamLAND) started operation on January 
21st 2002. The observation of i7 e 's from nuclear power reactors was the primary object of the 
KamLAND experiment because it was expected to solve the solar neutrino problem and set a 
strong constraint on the neutrino oscillation parameters. The search for solar anti-neutrinos was 
undertaken to examine two types of theoretical framework, spin-flavor precession and neutrino 
decay. This letter describes the reactor i7 e observation^ and the solar i7 e search™ KamLAND 
may also observe other neutrino sources. The observation of geo- neutrinos in high statistics will 
reveal quantitative information of U and Th in the earth. Observation of 7 Be solar neutrinos will 
be the next goal of KamLAND, however we still need to purify the liquid scintilator to achieve 
a detection of the Be neutrino signal. Other possible neutrino sources, such as supernova 
neutrinos and relic neutrinos are also an exciting prospect. KamLAND also has the ability to 
detect nucleon decay. 



1.2 Detector Overview 

KamLAND is located 1000 m under ground in the Kamioka mine in Japan. The detector is 
composed of a sphere shaped inner detector (ID) and a cylindrical shaped outer detector (OD). 
The inner detector which consists of 1200 m 3 of liquid scintillator (LS) (90% dodecane, 20% 
pseudocumene and 1.52 g/1 of PPO) is contained within a 13 m diameter balloon (134 /im thick 
transparent nylon/EVOH composite film) which sits inside 1800 m 3 of buffer oil (50% dodecane, 
50% isoparaffm). On the inside edge of ID sphere 1325 fast timing 17 inch diameter PMTs and 
554 20 inch diameter PMTs collect photons from the LS. Only the 17 inch PMTs were used for 



this analysis. The photo coverage was about 22%. The OD is filled with purified water and is 
instrumented with 225 20 inch diameter PMTs. The OD is a water cherenkov detector used to 
eliminate cosmic ray muons. It also acts as a shield against radioactivity from the surrounding 
rock. KamLAND is now stable and taking data 24 hours/days except for calibration runs or 
unexpected shutdowns. The average trigger rate is ~ 30Hz with the primary trigger threshold 
(200 PMT hits, ~0.7MeV). 

1.3 Energy and Vertex Calibration 

To calibrate the energy scale from lMeV to several MeV, various radioactive sources are de- 
ployed along vertical-axis of detector. The sources are 68 Ge (1.012MeV 7 + 7), 60 Co (2.506MeV 
7 + 7), 65 Zn (1.116MeV 7) and 241 Am/Be (2.20, 4.40, 7.6MeV 7). We also use 7-rays from neu- 
tron capture on proton (2.2MeV) and neutron capture on 12 C (4.95MeV). These neutrons are 
generated by spallation from the passage of cosmic-ray muons in the LS. The estimated energy 
resolution is about 7.5%/ y/ E(MeV). Light yield is 300 p.e./MeV. Radioactive sources are also 
used to obtain the positioning bias of the vertex reconstruction. The bias of the reconstructed 
vertex position was less than ±5cm along z-axis within the fiducial volume. 

1.4 Detection of i7 e by Delayed Coincidence 

Electron anti- neutrino is detected using the inverse (3 decay reaction in the LS, 

v e + V ~ > e+ + n - (1) 

The generated positron is immediately annihilated to 2 7 generating a prompt signal. On the 
other hand the neutron is thermalized in the LS and captured on a proton after ~200 (is, 
generating a 2.2 MeV 7 as a delayed signal. By requiring the coincidence of prompt and delayed 
signal, background can be reduced dramatically. The visible energy of prompt signal is related 
to the i7 e energy by: 

^vis = E v - e -{Am np + m e )-T n (9) + 2m e (2) 
= E v - e - 0.782MeV - T n {9) (3) 

where E„ e is energy of the z7 e , Am np is the mass difference between the neutron and proton, m e 
is the electron mass and T n (9) is the kinetic energy of the neutron scattered by an angle 9. The 
detection efficiency for inverse (3 decay events is estimated from Monte Carlo simulation and 
calibration data to be 84.2 ± 1.5%. The main contributions to the detection inefficiency are the 
cuts on the distance between the prompt and delayed vertices (89.8 ± 1.6%), the time between 
the prompt and delayed vertices (95.3 ± 0.3%), neutron capture on protons (99.5%), and the 
energy of the delayed event (98.9 ± 0.1%). 

1.5 Spallation Events after Muon 

Although KamLAND is located deep underground (2700 m water equivalent), muons originating 
from cosmic-rays pass through the detector at a frequency of 0.34Hz. An energetic muon can 
destroy a carbon nucleus in LS by spallation. Various radioactive isotopes and neutrons are 
generated by muon spallation. Tabled shows a list of generated isotopes. 

Spallation events due to muons are an important source of background. Most of the long 
life spallation products emit a single /3-ray so they are eliminated by the delayed coincidence 
requirement. However 9 Li and 8 He produce both a j3~ and neutron. For these we applied the 
following criteria to cut the correlated spallation events. (1) 2 ms veto is applied after any muon 
events. (2) Additional 2 s veto is applied for energetic muon events when the ionization energy 



Table 1: Radioactive isotopes generated by muon spallation in the liquid scintillator. 



Isotope 


T\/2 


#max (keV) 




20.2 ms 


13369 (/?") 


i2 N 


11.0 ms 


17338 (/?+) 


Xi Li 


8.5 ms 


20610 (/?") 


9 Li 


173.8 ms 


13606 (J3~, n) 


«He 


119.0 ms 


10653 (J3~, n) 


9 C 


126.6 ms 


16498 (/?+) 


"Li 


838.0 ms 


16006 (/?") 


e He 


906.7 ms 


3508 (J3~) 


«B 


770.0 ms 


17979 (/?") 



deposit is larger than 10 6 p.e. (~3 GeV). (3) For smaller energy deposits (less than 10 6 p.e.) a 
2 s veto is applied to events with vertices in a 3 m cylinder around the muon track. 

Muon spallation in the rock surrounding the KamLAND detector can generate fast neutrons 
that penetrate through the water of outer detector. These fast neutrons can cause a correlated 
background since they can generate both prompt and delayed signal in the LS via a recoil 
proton and neutron capture 7-ray. Adopting a 5m radius fiducial volume eliminates most of this 
background. 

Muon spallation provides not only background but also important calibration sources. Neu- 
tron captured on proton (2.2MeV 7), neutron captured on 12 C (4.9MeV 7) and 12 B decayed (5 
are used for energy and vertex calibrations. 

2 Result of Reactor v e Observation 

KamLAND is the first experiment to observe v e disappearance from a reactor source. - In this 
section the reactor experiment analysis is described. The dominant background of reactor z7 e 's 
is geo-i^ e 's from (5 decay of U and Th in the earth. The geo-i7 e flux has a expected energy 
distribution from 1 MeV to 2.4 MeV.^ Since the total flux of geo-z7 e 's nas large uncertainty, we 
applied an analysis energy threshold at 2.6 MeV for the reactor analysis. 

2. 1 v e Flux from Reactors 

Nuclear reactor power plants are strong electron anti- neutrino source. i7 e 's are generated by (3 
decay of daughter nuclei from the fission of various fuel components, such as 235 U, 238 U, 239 Pu 
and 241 Pu. The flux of v e from a reactor is calculated from the thermal power generation and the 
distance between the detector and reactor. The thermal power information is provided by the 
power companies with better than 2% accuracy. About 70GW (7% of the world total) of reactor 
power is generated at a distance of 175±35km from Kamioka. This corresponds to 80% of v e 
flux at KamLAND. Fortunately this distance is suitable to observe LMA neutrino oscillation. In 
our data set of 145.1 live days, the expected number of reactor z7 e events in the fiducial volume 
was 86 ± 5.6 events. The systematic uncertainties for reactor z7 e detection are listed in the table 
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2.2 Event Selection for Reactor v e 

The selection criteria for v e events is (1) total charge cut, less than 10,000 p.e. (~30MeV) 
with no OD veto signal and muon spallation cut. (2) fiducial volume cut (R < 5m), (3) time 



Table 2: Systematic uncertainties of reactor j/ e (%). 



Total LS mass 2.1 
Fiducial mass ratio 4.1 
Energy threshold 2.1 
Efficiency of cuts 2.1 
Live time 0.07 


Reactor power 2.0 
Fuel composition 1.0 
Time lag of (3 decay 0.28 
v spectra 2.5 
Cross section 0.2 


Total systematic error 


6.4 



correlation cut (0.5/xs < AT < 660/is), (4) vertex correlation cut (AR <1.6m), (5) delayed 
energy window cut (1.8MeV < E^ e i ay < 2.6MeV), and (6) cut on the delayed vertex position 
more than 1.2m from central vertical axis to eliminate background from thermometers of LS. 

After applying our cuts and a 2.6MeV analysis threshold, 54 events remained. The total 
background was estimated to be 0.95±0.99 events, where accidental background is 0.0086 ± 
0.0005, background from 9 Li/ 8 He is 0.94 ± 0.85 and fast neutron is less than 0.5 events. The 
ratio of the number of observed reactor i7 e events to expected events without oscillation is 

N bs ~ N BG = Q 6n ± .085(stat) ± 0.041(syst). (4) 

^expected 

This result indicates i7 e disappearance with 99.95% C.L. 
2.3 Interpretation with Neutrino Oscillation 

Neutrino oscillation is one of the most probable explanations to understand the observed deficit 
of z7 e 's. A "Rate analysis" was performed where oscillation parameters are examined by defining 
a 

2 (0.611 - R(sm 2 26,Am 2 )) 2 
Xrate ~ 0.085 2 + 0.041 2 d 

where, i?(sin 2 29, Am 2 ) is the ratio of expected number of events with oscillation to expected 
number of events without oscillation. Also a "Shape analysis", comparing the normalized energy 
spectrum of i7 e was performed. Figure ^ shows the visible energy spectrum of the prompt signal 
that corresponds to the i7 e spectrum. The upper figure is the expected reactor v e energy spectrum 
with contributions from geo-i7 e (model Ia)^ and accidental background. The lower figure shows 
the energy spectrum of the observed prompt events (solid circles with error bars), along with the 
expected no oscillation spectrum (upper histogram, with geo-j7 e and accidentals shown) and best 
fit (lower histogram) including neutrino oscillations. The shaded band indicates the systematic 
error in the best-fit spectrum. The vertical dashed line corresponds to the analysis threshold at 
2.6MeV. 

A maximum likelihood function was used to define the combined x 2 °f the rate and shape 
analysis. 

Xrate+shape = Xrate + Xbg( N BGI~2) + Xdist («1~4) ~ 2 log L shape (sm 2 6, Am 2 , N BG1 ^ 2 , <*1~4) (6) 

where L s h ape is the likelihood function of the spectrum including deformations from various 
parameters. Nbgi~2 are the estimated number of 9 Li and 8 He backgrounds and evi^4 are the 
parameters for the shape deformation coming from the energy scale, resolution, v e spectrum and 
fiducial volume. 

Figure |2] shows the neutrino oscillation parameter region for two neutrino mixing at the 
95% C.L. The excluded region from ChoozEl and Palo Verde^ are shown at the top region. 
The allowed region of the LMA solution from solar neutrino experiments^ is shown in the 



middle region. The KamLAND "Rate" analysis excluded most of the region except LMA. The 
allowed region of LMA is split into two regions by the KamLAND "Rate+Shape" analysis. The 
remaining lower region is called LMA-1 and upper one is called LMA-2. 
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Figure 1: Expected z/ e spectrum from the reactor (upper) and visible energy spectrum of the prompt signal (lower) 
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Figure 2: Neutrino oscillation parameters excluded by Chooz, Palo Verde, KamLAND (Rate) and allowed by 
solar-i/ combined (LMA), KamLAND (Rate+Shape). Two regions (LMA-1, LMA-2) remain at 9596.C.L. 



2.4 Future Prospects of the KamLAND Reactor u e Observation 

KamLAND is taking data continuously and an update on the reactor v e analysis will be coming 
soon. The expected v e flux is changing because some reactor power plants have shutdown for 
maintenance from September 2002. Figure shows the time variation of the expected i> e flux 
from reactors and the observed number of v e events at KamLAND. It is expected that a flux 
modulation analysis will check the consistency of the reactor z7 e deficit. From 2006 another 
strong reactor "Shika2" located 88km from KamLAND will begin operation. Figure 0] shows 
the contribution to the v e spectrum from "Shika2" considering the neutrino oscillation solutions 



LMA-1 and LMA-2. Since the distance 88km is the most sensitive region to distinguish LMA-1 
and LMA-2, it is expected that spectrum shape analysis will exclude one of these two solutions. 
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Figure 3: Time variation of the expected reactor v e flux at KamLAND is shown (solid line). The observed number 
of reactor v e events is also shown (solid circles with error bars). 
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Figure 4: Expected u e spectrum from the new reactor "Shika2" after 3 years exposure, considering the neutrino 

oscillation solutions LMA-1 and LMA-2. 



3 Result of Solar v e Search 



Although neutrino oscillation is the most favored solution to explain the sol ar n eutrino deficit 
other possible solutions are not completely excluded. The search for solar v^SlII j s an effective 
test to examine other exotic solutions. In this section an analysis of 185.5 live days KamLAND 
data is discussed. 

3.1 Possible Mechanisms to Allow the Solar v e 

It is generally believed that neutrinos from the Sun are generated by fusion of light nuclei (mainly 
protons) in the core. Anti neutrinos are not generated among the many processes of the fusion 
in the Sun, so we should assume some other mechanism to produce solar v e . In this letter two 
models are discussed. The first is a hybrid model^ of resident spin flavor precession (RSFP) 
and the Mikheyev-Smirnov-Wolfenstein (MSW) effect. In this model v e with non-zero transition 
magnetic moment can evolve into or y T while propagating through intense magnetic fields in 
the solar core and they can in turn evolve into i7 e via the MSW effect. The other mechanism 



comes from a model of neutrino decay^, where a heavy neutrino mass eigenstate may decay into 
a lighter anti-neutrino mass eigenstate. 



3.2 Event Selection for Solar v e 



The dominant component of solar v e flux above the current KamLAND analysis threshold 
(2.6MeV) is the 8 B neutrino flux, which extends up to 14 MeV. The reactor i7 e 's, who's en- 
ergies extend up to 7 MeV, become a source of background events for the solar v e analysis. 
Figure El shows expected visible energy spectrum of the reactor z7 e together with the v e spectrum 
from 8 B neutrino assuming that 1% of flux is converted from v e to i> e . The lower analysis thresh- 
old (7.5MeV) corresponds to the end point of reactor z7 e spectrum and upper one (14.0MeV) 
corresponds to that of the 8 B v e spectrum. The criteria to select solar v e events is (1) total 
charge cut, less than 10,000 p.e. (~30MeV) with no OD veto signal and muon spallation cut. 
(2) fiducial volume cut (R < 5.5m), (3) time correlation cut (0.5/is < AT < 660/xs), (4) vertex 
correlation cut (AR <1.6m), (5) delayed energy window cut (1.8MeV < Ed e i ay < 2.6MeV), 
(6) LS thermometer cut, and finally (7) prompt energy cut for solar i7 e (7.5MeV < E prompt < 
14.0MeV). Data from 185.5 live days were included in this analysis. Figure shows the prompt 
and delayed energy distribution of the candidate events before cut (7). After cut (7) no events 
remained. The total expected background was estimated to be 1.1±0.4 events, composing the 
backgrounds from reactor i7 e (0.2±0.2), atmospheric v (0.001), fast neutrons (0.3±0.2), acci- 
dental coincidences (0.02) and 8 He/ 9 Li (0.6±0.2). Systematic uncertainties are summarized in 
Table El 
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Figure 5: Expected visible energy spectrum of reactor v e together with the is e spectrum from 8 B neutrino assuming 
that 1% of flux is converted from v e to i7 e . The analysis energy range of solar z7 e was 7.5MeV< E pr0 mpt <14.0MeV. 



Table 3: Systematic uncertainties of solar v a (%). 



Detection Efficiency 1.6 
Number of target protons 4.3 
Live time 0.07 


Cross section 0.2 
Energy threshold 4.3 


Total systematic error 


6.3 



3.3 Upper Limit of Solar v e Flux 

Since no candidates are found in the corresponding; energy region, the upper limit of solar z7 e flux 
is calculated using the Feldman-Cousins method'^ with u e cross section (a = 6.88 x 10~ 42 cm 2 ), 
detection efficiency (e = 0.841), live time (T = 1.60 x 10 s) and number of target protons 
(p p x f v = 4.61 x 10 31 ). At the 90% C.L. the upper limit of solar i7 e flux is fo B < 3.7 x 10 2 
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Figure 6: Prompt and delayed energy distribution from the solar j/ e candidates before the final energy cut (7) 
(see text). No events remained in the energy range 7.5MeV< E prompt <14.0MeV and 1.8MeV< Edeiay <2.6MeV. 



cm 2 s . Since 29.5% of total 8 B neutrino flux of 5.05^Q[gJ x 10 6 cm 2 s 1 is contained within the 
energy window of this analysis, the v e to i> e conversion probability without neutrino oscillation is 
2.8 x 10~ 4 at 90% C.L. If we assume this co nver sion occurs by RSFP+MSW^ and also assume 
the recent best fit of oscillation parameters (sin 2 6 = 0.28, Am 2 = 7.2 x 10~ 5 eV 2 ), and a 
solar magnetic field model^l the upper limit of the neutrino transition magnetic moment /U„ 
and magnetic field B max is estimated as \i v ■ B max < 1.4 x 10~ 5 /ieG (90%C.L.). If we assume 
neutrino decajf^we can constrain the lifetime limit to Tijmi > 6.7 x 10 -2 s/eV. 

4 Summary 

KamLAND has observed an evidence for the reactor i7 e disappearance at 99.95% C.L. Assum- 
ing CPT invariance only the LMA solution is compatible with the deficit. We got an upper 
limit of the solar u e to v e conversion probability 2.8 x 10" 4 at 90% C.L. in the energy range 
8.3MeV < Ejj e 14.8MeV. The KamLAND experiment is supported by the COE program of the 
Japanese Ministry of Education, Culture, Sports, Science, and Technology and the United States 
Department of Energy. 
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